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The formation of a strange or hybrid star from a neutron star progenitor is believed to occur when 
the central stellar density exceeds a critical value. If the transition from hadron to quark matter is 
of first order, the event has to release a huge amount of energy in a very short time, and we would be 
able to observe the phenomenon even if it is at cosmological distance far from us. Most likely, such 
violent quark deconfinement would be associated with at least a fraction of the observed gamma 
ray bursts. If we allow for temporal variations of fundamental constants like Aqcd or Gn, we can 
expect that neutron stars with an initial central density just below the critical value can enter into 
the region where strange or hybrid stars are the true ground state. From the observed rate of long 
gamma ray bursts, we are able to deduce the constraint Gn/Gn ^ 10~^^ yr~^, which is about 5 
orders of magnitude more stringent than the strongest previous bounds on a possible increasing Gn ■ 



Introduction - Search for temporal variation of 
"fundamental constants" has attracted a lot of interest 
in last years [ij. Such a phenomenon is quite a gen- 
eral prediction of many different frameworks (superstring 
theories, scalar-tensor theories, models in extra dimen- 
sions, etc.) when the dynamics of the whole universe is 
taken into account and a clear sign of violation of the 
Strong Equivalence Principle (SEP) [2]. On the other 
hand, the Einstein Equivalence Principle (EEP) - which 
is the essence of the geometrical theory of gravitation - 
refers only to non-gravitational physics and thus allows 
the Newton constant Gn to be time dependent, but still 
forbids any temporal variation of non-gravitational pa- 
rameters (a, Aqcd, Gp, etc.) Q. 

Since quantitative predictions are impossible at 
present, the common approach is based on model inde- 
pendent investigations of possible temporal evolution of 
physical quantities that should instead be constant in 
the standard theory and constraints are often reported 
assuming variations that are linear in the cosmological 
time: indeed one can expect to expand any false con- 
stant in a power series of time and that the linear term 
is the leading-order correction. For this reason, essen- 
tially all the works in the literature on this subject do 
not specify any theoretical framework and the discussion 
is at a pure phenomenological level. However, this ap- 
proach cannot be always justified and hence, in general, 
one cannot directly compare bounds obtained today in 
laboratory with others coming from arguments involving 
the physics of the early universe. 

The topic of time varying "fundamental constants" has 
received additional interest after the claims of some ev- 
idence of time variation of the fine structure constant 
a [3| and of the proton to electron mass ratio nip /me v^. 
There are indeed good theoretical arguments to believe 
that a possible temporal evolution of a implies time shift 
for other quantities as well Q . However the observational 
situation is quite controversial, since other groups have 
not confirmed the results Q. 

A relevant quantity which sets strong interactions and 



hadron masses is the QCD scale Aqcd- The most strin- 
gent constraint on its possible temporal variation comes 
from an analysis of isotopic abundances of the Oklo ura- 
nium mine in Gabon, a natural fission reactor that oper- 
ated about 2-10^ years ago for about 10^ years. Assuming 
a hncar time dependence, one gets [3| 



Aqcd/Aqcd 



< 10 



'17 



yr 



(1) 



Other constraints on Aqcd are usually quite model de- 
pendent and, in any case, at least 2 orders of magnitude 
weaker than the bound in P]). 

However, among all the "fundamental constants" a 
privileged role is played by the gravitational constant 
Gn, whose variation violates only the SEP. The tight- 
est constraint in the literature on a linear temporal evo- 
lution of Gn is provided by the Lunar Laser Ranging 
experiment, which has monitored Earth-Moon distance 
for about 30 years. The result is Q 



Giv/G'jv = (4±9)-10-^^yr 



(2) 



On the other hand, the earliest event in the history of 
the universe which has left a reliable record of the past 
value of Gn is the Big Bang Nucleosynthesis: from the 
primordial abundance of "^He we can conclude that, when 
the universe was only 1 s old, i.e. ~ 13 Gyr ago, Gn could 
not differ more than about 10% from its present value [^. 
Other interesting constraints on Gn/Gn, at the level of 
lO^^'^yr^^ in the case of linear time evolution, can be 
deduced for instance from the orbital period of binary 
pulsars 10], from the comparison of masses of old and 
young neutron stars [11] and from the so called gravito- 
chemical heating in neutron stars il^. 

In this letter, we consider the effects of a temporal vari- 
ation of Aqcd and Gat on the equilibrium configuration 
of compact stars: indeed, non-zero Aqcd or Gn could 
change the true ground state of some stars and induce 
conversion of neutron stars into strange or hybrid ones. 
Here the keypoint is that the transformation probably 



involves a first order phase transition, leading to the re- 
lease of a large amount of energy in a short time. This 
makes even small variation of Aqcd or Gn to produce 
spectacular events, observable at cosmological distances. 
Taking into account the possible events which could be 
somehow associated with these phenomena, we can esti- 
mate an upper limit for the decreasing rate of Aqcd or 
for the increasing rate of Gn- 

To avoid any possible misunderstanding, a comment is 
in order here; we can only constrain dimensionless quan- 
tities, such as a. When we talk about temporal evolu- 
tion of dimensionful constants, like Aqcd and G^, we 
implicitly assume a particular system of units and again 
we constrain the dimensionless ratio to standards chosen 
as units and constant by definition 'l|. This is partic- 
ularly evident for example in some scalar-tensor theo- 
ries of gravity, where one can work in the Einstein frame 
(where Gn is constant) or in the Jordan frame (where 
particle masses are constant and Gn can change), see 
e.g. Ref. [1J|. However, this is a more general result and 
one can always choose an arbitrary dimensionful quan- 
tity as standard unit and then compare it with other 
quantities. And this is basically also what we do here: 
when we assume Gn constant, we allow for variation of 
Aqcd, and when we assume Aqcd constant we discuss 
the possibility of non-zero Gn-, but we compare always 
the QCD interaction with the gravitational one. For this 
reason, the possibility of a compensation between the 
temporal variation of the QCD energy scale Aqcd and 

— 1/2 

of the gravitational energy scale Mpi = Gjy has no 
physical meaning: only their dimensionless ratio is an 
observable. Of course, here we cannot exclude the pos- 
sibility that both Aqcd and Mpi change exactly in the 
same way with respect to another dimensionful constant, 
for example the vacuum expectation value of the Higgs 
field V, but such a case sounds quite ad hoc and would 
require an unwanted fine-tuning. 

Compact stars - Compact stars are the end-product 
of heavy stars after supernova explosion [ij, llH . Their 
mass is typically ^ 1.5 M© and their radius is about 10 
km, so that they are extremely dense objects: in partic- 
ular we expect that the central baryon density, pc, is at 
the level of 2 - 10 times normal nuclear matter density 
po = 0.16 fm~'^. At such high densities, matter would 
be made of neutrons, with a smaller fraction of protons 
and electrons. However, in the upper region of this inter- 
val, the ground state is probably represented by strange 
matter of m, d and s deconfined quarks: the latter are 
fermions and introducing a third flavor there are new low 
energy available states which can reduce the total energy 
of the system. If so, compact stars with a central den- 
sity higher than a critical value pcrit would be strange 
stars (stars made entirely of strange matter) or hybrid 
stars (stars with a strange matter core and an outer part 
of hadron matter), depending on QCD physics and/or 
their mass. For a review, see e.g. Ref. [15|. Thanks to 
these extreme conditions, impossible to reach in labora- 
tories on the Earth, compact stars can really be seen as 



a unique opportunity to investigate new physics [Ig . Il7l | . 
The results of this letter are essentially based on the 
following two assumptions: 

1. Strange or hybrid stars exist, but not all the com- 
pact stars are strange or hybrid. 

2. The transition from hadron to strange matter is of 
first order. 

The first statement means that we assume that the pop- 
ulation of compact stars is made of both neutron and 
strange/hybrid stars. In other words, the critical density 
Pcrit has to be in the range of the typical central density 
of compact stars, so that the ones with pc < Pcrit are 
neutron stars, while if pc > Pcrit they are strange or hy- 
brid. As for the second hypothesis, it is indeed what is 
commonly believed, i.e. at high baryon density the phase 
transition from hadron to strange matter seems to be of 
first order 18]. Unfortunately, up to now no model inde- 
pendent argument supporting that hypothesis has been 
provided. 

Let us now introduce the information which is needed 
to constrain the temporal variation of the QCD scale and 
of the gravitational constant. From considerations on 
stellar evolution, we can say that about 0.1% of the stars 
in galaxies are compact ones, that is neutron or strange 
stars (see e.g. Ref. [14|). This implies that there are 
about 10^ compact stars in a typical galaxy and that 
the whole visible universe contains a total number of 
Ntot ~ 10^° [13] • As for their central baryon density 
Pc we assume for simplicity that compact stars are uni- 
formly distributed with pc in the range (2 p^, 10 po), that 
is the central density distribution is 



n{Pc) 



Nt, 



>Po 



(3) 



for 2 Po < Pc < 10 Po and otherwise. Here we would 
like to stress that this simple choice is not crucial for our 
conclusions. Indeed, we could consider another distri- 
bution and get roughly the same result if the probabil- 
ity of finding a star with a given central density is the 
same within an order of magnitude: for example, Eq. ([3]) 
could be seen as the approximation of a gaussian of height 
~ Ntot/S>po and width r^ 8po. Following our first work- 
ing hypothesis, 2 po < Pcrit < 10 po and the possibility 
of temporal variation of Aqcd and Gn may induce the 
transformation of neutron stars into strange or hybrid 
ones, because pcrit depends on Aqcd and pc on Gn- In 
addition to this, since we take the phase transition to be 
of first order, in the conversion the compact object has 
to release a huge amount of energy in a very short time, 
say 10 s. This energy is the difference of binding energy 
between strange/hybrid star and neutron one. All the 
models agree with a typical energy release at the level of 
lO^'^ — 10^'^ erg [i^. This huge amount of energy can be 
converted, at least in part, to photons and e~^e~ pairs. 
Various mechanisms have been proposed [20| , all of them 



producing a gamma emission with an energy easily ex- 
ceeding 10^° — lO'^^ erg. Such a huge power in 7 can be 
at the origin of at least a fraction of the so called long 
Gamma Ray Bursts (GRBs) [2l|. 

Temporal variation of A-qcd - A temporal vari- 
ation of Aqcd must induce a change in the value of the 
critical baryon density pcrit- A rough estimate can be 
obtained as follows. We start by considering two types 
of quar k models, the bag-like [24I and the NJL-like mod- 
els [23] ■ In bag models, the crucial parameter regulating 
Pcrit is the pressure of the vacuum B, whose numerical 
value is independent of the presence of the quarks. There- 
fore the value of B is controlled by the only dimensional 
parameter present also in pure gauge QCD, i.e. Aqcd- 
One can therefore write 



Pc 



u^aK^ 



QCD 



(4) 



where a is a dimensionless order one coefficient. In chiral 
models, and in particular in the NJL-like ones, gluons 
are not present as dynamical degrees of freedom. It is 
anyway possible to introduce a quantity whose physical 
meaning is similar to the one of B in bag models (see 
e.g. Ref. [231). The numerical value of that quantity 
is dominated by the value of the constituent masses of 
the quarks. In NJL-like models the mass of the strange 
quark is so large that it plays a minor role compared to 
the up and down quarks. Since the constituent masses of 
the u and d quarks are almost independent of the small 
value of the current mass, one can conclude that also 
in chiral models the "pressure of the vacuum" depends 
almost only on Aqcd- 

In the bag model, the relation between pcrit and B has 
been evaluated in a large number of calculations and is 
typically of the order of 15, 24] 



0.003 ^ „ 

^ MeV 



(5) 



where Apcrit and AB are respectively a variation of 
Pent and B. From Eq. (g]) we find that AB/B = 
4AAqcd/Aqcd and hence 



Apc 



IOPqAAqcd/Aqcd- 



(6) 



Let us note that A Aqcd < implies Apcrit < 0: a 
time decreasing Aqcd makes pcrit decrease as well and 
neutron stars with a central density just below the crit- 
ical value may convert into strange or hybrid ones. The 
rate of strange/hybrid star production due to the possible 
temporal evolution of Aqcd we can observe is 



N,. 



-en{pc) Pcrit 



10 



20 



A 



QCD 



A 



QCD 



(7) 



Here e is the fraction of compact stars in the visible uni- 
verse whose conversion from hadron to strange matter 
can be seen from the Earth and the equation is valid 
only for pc < (hence Ngg > 0). If Pc > 0, strange 
stars should convert into neutron ones, but in this case 



we do not expect any violent phenomenon: more prob- 
ably the picture would look like a slow evaporation of 
strange matter into hadrons. Because of this, no bound 
can be deduced for the case Aqcd > 0. 

In conclusion, according to our theoretical framework, 
a time decreasing Aqcd must be accompanied by spec- 
tacular events, that is the formation of strange or hybrid 
stars from neutron star progenitors. As already men- 
tioned, the only known events which can be associated 
to such violent phenomena are the GRBs [2l| . In partic- 
ular we are interested in the so-called long GRBs, which 
are known to be associated with the gravitational col- 
lapse of massive stars. Here we would like to stress that 
we are not assuming that all or some of the long GRBs 
are related to the formation of a strange star, but we are 
saying that the conversion from a neutron to a strange 
star, if it occurs in the universe, is certainly a quite par- 
ticular event which likely releases a huge amount of en- 
ergy (~ 10^^ erg) in a short time and that at present we 
have not yet observed other phenomena compatible with 
this picture but long GRBs. From this point of view, we 
can assert that the rate of strange star production can- 
not exceed the one of GRBs. Since today the observed 
GRBs have a mean redshift z ~ 3 and are more or less 
uniformly distributed up to redshift z « 5, we would be 
able to observe most of the events in the whole visible 
universe associated with the transition from a neutron 
star to a strange star. This means that we can put e '-~^ 1 
in Eq. ^. In the end, since the observed rate of GRBs 
is 



N, 



GRBs 



3,,^-l 



lO'^yr 



(8) 



we deduce the following constraint on a possible linearly 
time decreasing QCD scale [28 1 



Aqcd/Aqcd<10 ''yr 



(9) 



One caveat is in order. If the transition from a neu- 
tron star to a hybrid or a quark star produces a beamed 
emission, the actual limit on Ncrbs is larger. GRBs 
are known to be beamed, and the correction to the 
rate due to beaming is of the order of 75 ± 25 [25J . 
If the same beaming factor applies to GRBs associated 
with quark deconfinement, than Nqj^bs ~ lO^yr"^ and 
-Aqci?/Aqcd< 10-15 yr-l. 

Temporal variation of Gn - Let us now consider 
the possibility of a non-zero Gn- Now the critical den- 
sity Pcrit is unchanged, because it is set by QCD physics. 
However, like all the other stars, neutron stars are self- 
gravitating system whose configuration is determined by 
the balance between the attractive gravitational force 
and the particle pressure. Hence a variation of Gn in- 
duces a change of the central stellar density. Of course, 
if a neutron star has a central density just below the crit- 
ical value, an increase of Gn makes pc pass the critical 
value and the neutron star can convert into a strange 
one. A rough estimate for the bound on Gn /Gn can 
be deduced from the following simple considerations. As 



stellar model, we can take a system of N non-relativistic 
ferm.ions of mass m which interact only gravitationally. 
The energy per unit volume is 



e ~ - 



3p% 



5 2m 



mp--- 



SGNm^P^V 



R 



(10) 



where p is the particle number density (for the sake of 
simplicity we assume it is constant in the star), pp ~ 
{Stt^pY^^ is the Fermi momentum, R the stellar radius 
and V the stellar volume. The first term on the right 
hand side of Eq. pU]) is the particle kinetic energy per 
unit volume, the second term is the particle rest energy 
per unit volume and the last one is the gravitational po- 
tential energy per unit volume. If we multiply Eq. pO|) 
by 1/yO and then we replace p with N/V, we get the en- 
ergy per particle as a function of the stellar radius R, 
since N must be constant. It is straightforward to find 
that the equilibrium radius is proportional to 1/Gn. In 
the Newtonian framework, for a star with constant den- 
sity we have p oc R~^ oc G^ . At first approximation, we 
can expect the same dependence of the central density of 
compact stars on the gravitational constant Gn and the 
relation between a small change of Gn and pc is 



Apjp, «3AGjv/G 



N 



(11) 



Hence, replacing pcrit with —pc in Eq. ([7|), we find the 
bound on the temporal evolution of Gjv 



Gjv/GAr<10-"yr-i. 



(12) 



As we have already said for the case of the QCD scale 
Aqcdj such a bound is applicable only in one direction, 
here for an increasing gravitational constant. However, 
Eq. (I12p is about 5 order of magnitude more stringent 
than the present best constraint of Eq. ([2]) . Also in this 
case, if the gamma emission is beamed the limit becomes 
less stringent, Gn/Gn S lO^^'^yr^^ 

Conclusion - In this letter we have discussed the 
possibility that the QCD scale Aqcd is decreasing in 
time or that the Newton gravitational constant Gn is 
increasing, considering their effects on the equilibrium 
configuration of compact stars. Under a set of reasonable 
assumptions, we have found the constraints 



The possibility of time varying "fundamental con- 
stants" is not exotic, but instead quite a general pre- 
diction of most theories of gravity beyond general rel- 
ativity. Assuming that compact stars can be neutron 
or strange/hybrid stars, depending on the value of their 
central density, we have to expect that Aqcd < or 
Gat > induces the transformation of some neutron stars 
into strange or hybrid ones. Then, if the phase transition 
from hadron to quark matter at such high density is of 
first order, as it is commonly believed, the energy release 
is huge and we should be able to observe the event even 
if it is at cosmological distant far from us. At present, 
the only phenomena compatible with this picture are the 
GRBs and, from their rate, we arrive at the constraints 
on the temporal evolution of Aqcd and Gn- On the 
other hand, if the origin of GRBs were completely differ- 
ent, it would mean that at present we do not know any 
phenomenon which can be associated to the formation of 
a strange/hybrid star from a neutron star and our bounds 
would become much stronger. Of course, it is important 
that future investigations will confirm our working hy- 
pothesis. In the end, we would like to stress that our 
constraints Q and P^ (or p^ and p^ ) are very com- 
petitive, even from the point of view of the covered time 
interval (at present, only BBN and CMBR arguments 
can probe earlier events in the history of the universe). 
In particular for the case of the gravitational constant, 
our bound on a linear time increasing Gn is much tighter 
than all the previous bounds, even if a possible beam- 
ing of the GRB associated with quark deconfinement is 
taken into account. On top of that, the previous bounds 
will be hardly significantly improved in a near future for 
the coming of new measurements and observations. On 
the contrary, our constraint could become much stronger 
once the origin of GRBs will be better known. 

Finally, it is tempting to discuss a recent analysis in- 
dicating an excess of long GRBs at large red-shift [2a |. 
More precisely, an evidence has been found of a GRB rate 
4 times larger than the one predicted assuming the GRB 
rate to follow the star formation rate. A possible inter- 
pretation of this result can be based on a non-linear evo- 
lution of the fundamental constants, with a more rapid 
variation in the young universe. 



Aqcd/ Aqcd < 10 

Gn/Gn 



-i^r-i 



< lO-^Vr^' 



(13) 
(14) 
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Indeed, the number of stars in a general galaxy can be 
estimated to be about 10^^, while the number of galaxies 
in the visible universe is roughly lO'^^. 
Here, as well as in Eq. (|12p below, the temporal variation 
is respect to the cosmological time, that is the time coor- 
dinate of the metric describing the expanding universe. It 
must be so because the prediction of time varying "funda- 
mental constants" arises when we consider cosmological 
solutions of theories beyond general relativity. 



